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The problem of the metastable phase decay has some specific solutions 
depended on the content of the condensating system. Namely, when the 
system contains only one type of heterogeneous centers one can successfully 
apply the iteration procedure 0. When there are at least two types of 
heterogeneous centers all attempts to apply the iteration procedure in 
the general situation fail. Then one has to consider the set of characteristic 
situations 0. This case will be the subject of our analysis. 

In some situations of the condensation process one can use the iteration 
procedure, but the reconsidered variant. In other situations one has to ap- 
ply another theoretical approaches 0. In it was announced that the 
mentioned characteristic situations cover all variants of the possible external 
conditions (and the content of condensated system). The concrete proof is 
very long and dull and it is omitted in 0] due to lack of volume. 

Ordinary the total number of heterogeneous centers is one and the same 
for all types of centers. For example, the total number of the positive ions 
is ordinary equal to the total number of the negative centers. This condition 
simplify the analysis and leads only to two characteristic situations. 

These situations are quite analogous to those considered in the general 
situation 0. Moreover, the way of proof is also quite analogous. That's why 
we restrict ourselves by the demonstration of the proof in this situation. 

To formulate two characteristic situations we have to recall some defi- 
nitions used in 0. We shall start with fundamental characteristics of the 
condensation process. We shall mark by lower indexes + and — two types of 



1 



heterogeneous centers. In the real situation of ions these indexes can mark the 
positive and negative centers correspondingly. The absence of index means 
that the value can be refered to both types of heterogeneous centers. 

To describe the first characteristic situation we denote by AF the height 
of the activation barrier. Certainly it is the function of the supersaturation 
( defined as 

where n is the molecules number density of the vapor and rioo is the molecules 
number density of the saturated vapor. As far as we have two types of centers 
we have two values of the activation barrier heights A+F and A_F. 

The first characteristic situation is the situation of the "strong unsym- 
metry" . It is characterized by the small value of parameter inverse to 

5AF= |A+F- A_F| 

More rigorously speaking we have to mention that really the small value 
of parameter 

exp(-|A+F- A_F|) 

is required. 

The second characteristic situation is the situation of the " moderate un- 
symmetry" . To formulate it we have to introduce the values of 

dC 

As far as there are two types of centers we have 

dA+F 



+ 



-c- 



and 

r_ — (- 

^ dC 

The second characteristic is characterized by the small value of parameter 

ir^-r_| 



Here we shall show that these two situations exhaust all possibilities. 
The structure of consideration is the following: 
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• In the first part we shall present general reasons which allow to hope 
that the mentioned situation can cover all possibilities of experimental 
conditions. 

• In the second part we shall derive the required overlapping for several 
models of heterogeneous centers. This illustrates the proof and gives 
the answer in some limit situations which will be used later. 

• In the third part we shall give the derivation of the statement for the 
case of ions. The way of derivation will be also spread to the general 
case. 

All definitions from |^ are acceptable. We shall follow the system of units 
and definitions used in these papers. 



1 General remarks 

The difference between two types of heterogeneous centers is induced by some 
abstract charge q. In the situation of ions it is a real electric charge, in other 
situations it is the abstract charge. 

Let us suppose that we can vary q starting from the zero valueQ. When 
g = we have A+F = A_F and r+ = r_. When q is increasing then 
and 5/S.F are increasing also. Let us stop when 



5AF = 1 (1) 



We mark this value of q as qh- 
Having introduced parameter 

6rAF = 



SAF 



A+F + A_F 
one can note that 



SrAFU « 1 (2) 



^In the case of ions there is an elementary charge of an electron but we use this way 
only for a formal derivation. 
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The activation barrier height AF can be presented as 



AF 



/•Co r 



where Co is the supersaturation when the activation barrier disappears. One 
can present the last integral as 

T(C') 

AF(C) = -p(Co-C) 
where (' is some value between ( and Co- Then 

c 

where 

T(C') C 

AC-(Co-C)^^ 

The value A^ has the sense of the characteristic distance from ( until the 
value where the essential activation barrier disappears in comparison with 
initial value (this value isn't Co )• The function (Co — C)/^C is a smooth 
function of C- 

Two types of heterogeneous centers induces two values A+C and A_C- 
One can easily show that 

A+C|,<,. ^ A_C| (3) 

We shall show the last estimate very qualitatively. 

Really, the barrier character of condensation implies that 

A+F> 1 

A_F> 1 

This leads also toPl 

^ rfA+F r+ ^ 

— — = — > 1 

^c c 



^It is necessary for continious description of the nearcritical region. 
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dA_F r_ 

Then the violence of the required condition leads to the violence of (|l|). 
So, one can now directly see (|^). 

One can keep in mind that the most sharp function of the supersaturation 
is r and the smooth functions are AF and A^. Then it is reasonable to 
transform this picture into the dependence on q. Namely, one can consider 
that for the given supersaturation ( the value F is the sharp function of q and 
the values AF and A( are more smooth functions. This can not be rigorously 
proven but seems to be a reliable qualitative picture. In the second section 
the similar facts will be justified for concrete types of heterogeneous centers. 

The mentioned approximate coincidence of A^( and A^( allows to in- 
troduce 

AC = ^(A+C + A_C) 

and approximately substitute A^^ and A_( by A^. 
Then one can rather approximately show that 

|A+F- A_F| ^ I— - — |AC (4) 



and as far as 



one can come to 



AC/C « 1 

|AC^^- ^AC 

where (p = (C+ + C_)/2. In the last relation one can take as approximately 
both and C_ with a rather small relative error. 

Now one can express A( through AF and substitute it into the last 
relation. Then it comes to 



± 



The last estimate solves the problem of overlapping of the mentioned 
regions. Really, as far as A^-F ^ 1 we see that at qh where | A+F — A_F| ~ 1 
the small value |r+ — r_|/r-|- is guaranteed. 
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We have to note that the vahdity of (0) is the matter of question. Cer- 
tainly, one can adopt A±F ~ r±^, but when we coming to the difference 
r_|_ — r_ the relative error increases many times. This disadvantage leads to 
some further remarks. 

2 Model systems 

Here we shall consider three simple models and show the necessary over- 
lapping directly. This will illustrate that the overlapping of the regions 
|A_|_F — A_F| > 1 and 6rT <^ 1 is rather natural. 

2.1 Pseudo homogeneous model 

Suppose that Rc + = Rc - = Rc q=o where R is the radius of the embryo, 
index "c" denotes the critical embryo. This corresponds to the weak influ- 
ence of the nuclei on the surface region of the nearcritical embryo. Later 
we shall approximately suppose that z/c + = z^c - = ^'c hom where u is the 
number of the molecules inside the embryo and index hom corresponds to 
the homogeneously formed embryo. 

Later the index hom differs from the subscript q = (Index "q=0" sup- 
poses only that the terms depended on the sign of the charge are put to zero. 
All other terms depended on the absolute sign of the charge are conserved.) 
The role of the terms depended on the sign of a charge is rather small in 
comparison with the role of the terms depended on the absolute value of a 
charge^. So, one can use the decomposition starting from g = when all 
terms independent on the charge sign taken directly into account. 

Now we shall return to direct calculations. For r+ and r_ we have 

r+ ~ l^hom — l^e + 

Then 

■^The homogeneous nucleation rate isn't between the rates of embryos formation on the 
"positive" and "negative" heterogeneous centers. Both positive and negative centers are 
the active centers of the condensation. 
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For |A+F — A_F| we have a very rough estimate 

|A+F - A_F| ~ F_(i/e -) - F+{ue +) 

All dependence on sign is in the last difference. We have to estimate the 
last difference by the smoothest dependence. Let us take the homogeneous 
dependence for this value. Then 

-) - F+(^e +) - ^{v'Jl - vT-) ^ P 

where a is the renormalized surface tension. 

Having estimated Uhom — ^e± > ki'e± with some parameter k close to 1 
one can use for (r_|_ — r_)/r-|- very rough (and smooth) estimate 



If S is small then 



where 



p = ^{u'JUl - {1 - kd f') ^ ASk 



4 = ^2 2/' ^ 1 



The property Ak ^ 1 guarantees that the regions (5 -C 1 and p > 1 are 
overlapped. 

2.2 Linear approximation model 

We shall start from the rigorous formula 

When C = (o± the difference in the r.h.s. goes to zero. Now we shall introduce 
approximation 

r± ~ 7±(C - Co±) 

which implies that the difference I'c ± — ± is the linear function of the 
supersaturation. This approximation has to be valid at some effective super- 
saturations which makes the main contribution into the activation barrier 
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height^ As it will be seen in the next subsection this approximation isn't 
valid when ( is near ^o- 

Having integrated the suggested approximation we come to 

For q = Qh we have 

i = Iy(C-Co+)^-^(C-Co-)^I 



or 



1 = A+F|l - ^( 



7-X - Co -x2| 



7+ C - Co + 

The value (C — Co -)/(C ^ Co +) has to be close to 1 or it has to be A+F — 
A_F ~ A+F ^ 1 which solves the situation. 
Then 

1 = A+F|l - —I 

7+ 

or 

1 > (A+F)-i = I 

7+ 



It can be valid only when 

11 - 

7+ 



1 - — I < 1 



One can see that 

7+ r± 



So 



,r+-r_, , 

1^ « 1 

which shows the overlapping. 

*Due to one can imagine AF as the result of growth of AF{() where C is 

falling from to 
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2.3 Moderate behavior model 



Now we shall construct the general model corresponding to the formation and 
disappearing of the metastable state. The simplest form of the free energy 
corresponding to the appearance of the gap in the region of the small sizes 
is following 

F ~ —b{u — vq) + c(z/ — 

where uq is the characteristic value|. 

Here b plays the role of the chemical potential (or supersaturation) and 
c is some negative parameter associated with the nuclei and independent 
(weakly dependent) on the supersaturation. 

Denoting z/ — z/q via x one can get 



AF 



4 63/2 



33/2 Id 1/2 

dAF 



2xn = Vr — Vp 



db 

which confirms F ~ z/c — z/g directly. 

The value of b is the variable, the value of c is supposed to be parameter. 
Now it is clear that r+ essentially differs from r_ only when c+ essentially 
differs from c_. But it means that A+F ~ |c|~^^^ essentially differs from 
A_F. As far as A±F > 1 we see that it means that |A+F - A_F| > 1. So 
the overlapping here can be also observed. 



3 Real systems 
3.1 Ions 

Now we shall investigate the case of ions. The free energy of the embryo 
formation on the electric charge q can be presented in leading terms as fol- 
lowing 

F = -bu + az/2/3 + cz/i/3 + (C2 + c,)z/"^/3 

^Here we are interested only in behavior of F near and the asymptotic behavior of 
F isn't essential (it is wrong). 
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where b is the excess of the chemical potential in a liquid phase, a is the 
renormalized surface tension, c and C2 are some coefficients depended on the 
absolute value of the nuclei charge, Cg is the coefficient proportional to the 
charge and, thus, depended on the sign of the charge. 

It is more convenient to use instead of u the variable p = v^l^ which leads 

to 

F = -bp^ + ap^ + cp+ (c2 + Cg)p~^ 
For the critical size one can get 

-36p^ + 2ap^ + cp^ = C2 + 

We shall present the critical characteristics in the following form 

i^c ± = i^c + 5±Vc 

l^e ± = t^eO + S±Ue 
Pc ± = Pc + S±Pc 
Pe ± = Pe + S±Pe 

where index marks the values when Cq — (but c and C2 are conserved). 
For po (both for critical and equilibrium values) we have 

-36po + 2apo + cpo = C2 

Then for 6p (here will be 6-^-p = —6p = 6p in the main order for both 
critical and equilibrium values) one can get 

-36(po + Spy + 2a(po + Spf + c(po + 5pf = C2 + 
and in the main order 

^ ~ -mpl + &apl + 2cpo 
For the critical value of the free energy one can get 

F ^ Foipo) + i^o (Po)^ + c,(po + Sp)-' 
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where 

Fq = -bp^ + ap^ + cp + C2p~^ 

F^' = -66p + 2a + 2c2^ 

Po 

Then one can come to 

+ - _ = 2c,Po 1 
+ - _ = 2c,Po e 

and 

A , F - A_F = 2c/ ^ ^ ■ 



Po c Po e 



Now we shall turn to get T = 1/^ — I'e- In the main order r+ = r_ = Fq 
z/co — i^eo- In the first order 



and 



r± = To ± (3p^5pc - 3pg(5pe) 



r+ - r_ 6p^5pc - 6pe5pe 



To Pc - Pe 

Then one can take for the last value the following estimate^ 

r+ - r_ 66pc 

r\-/ 

To Pc 

When q qh we get 

1 1 ^ 2c, 



1 = 2cg( 



and we see that Cq ^ 1. 

Then for + we can find 

To 



Po c Po e Po 



r+ - r_ 3po 



To Po c't^(Pco) 



^If we take here Spe instead Spc all consideration cen be repeated even in details. Then 
we have to take V = V{peo) ^ 1- 
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where function V is given by 

V{p) = -12bp^ + 6ap^ + 2cp 

As far as V{pco) ^ 1 we see that 

r^-r. 



To 

which proves the overlapping. 



< 1 



3.2 Generalization for the arbitrary system 

Now we can use the last constructions to investigate more general situation. 
Suppose we have F_,_, F_ and Fq. Then 

F± = Fo + 6±F 

For the critical and equilibrium values we have 

To find 5±fc,e one can use 



= 



or 



or 



which gives 



du du 

dFo . d ,dFQ.. . , db±F 



Ordinary in the leading term S±F — ±SF. Then 

dSF \ 

r L dv l'^='^e,e, 

()±I/c,e - ±^2^ 

(il/2 \l^=l^c,e, 
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Now we shall find A^F — A_F. For Fiur) we eet 



Then 

F . - F = 

3 " 



F^^-F^. = l^5vl + 25F\ 



and 



The first term corresponds to the opportunity missed in the previous section. 
Here we sgall consider it more correctly. 
When q = qh we get 



or 



1 (dSF\3 



To find r± we shall use T ^ — Ue- Then 



To l^c ~ l^e 

and 

- r_ 



To t'c ~ i^e 

the nearcritical region. Then 



For further analysis one can express in terms of the halfwidht AcZ/ of 



where AcZ/ is the halfwidth of the nearcritical region. Then (^) transforms 
into 

Id^Fo d6F 3 6 



3 (ii/^ c?!/ 



(A,z/)^ + 25F|.=.„^ + ...|.=.„^~l 
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One can easily prove that in the Lh.s. of last relation there is no com- 
pensation and get the estimates 



I I'D 



Then as far as 



3 du^ du 



J'(Ac^)' < 1 (6) 



To l^c - Ve 

according}] to ® one can see that 

r+-r_ 1 .3 

Having used the estimate 

du"^ du^ 

one can get 

|r+-r_| 



1/3 



To i^c - 

which solves the problem when the denominator isn't too small. But the 
small value of denominator is already investigated in the section "Moderate 
behavior model". 

For 6F ~ 1 we have 

I — aAF( ) 

To Vc - Ve Vc Ve 

which solves as far as here 5F ~ 1 the problem when z/^ isn't very close to v^. 
But this situation is already investigated in the section "Moderate behavior 
model" . 

The leading property which allows to justify all necessary estimates is 
the fundamental condition z/c ^ 1 which is necessary for the thermodynamic 
description of the embryo. 



^ If we choose another opportunity 5F ^ 1 then having adopted 5F ^ Cq/v " with 
some parameter a we repeat the previous section. 
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